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WAVENUMBER  -  FREQUENCY  STRUCTURE  OF  TURBULENT 
FLOW  IN  A  CHANNEL  FROM  A  DIRECT  SIMULATION 

1.  INTRODUCTION 


In  this  work  we  pre-sent  details  of  the  wavenumber  -  frequency  structure  of  tur¬ 
bulence  in  a  channel  obtained  from  a  direct  numerical  simulation.  The  work  parallels, 
in  certain  respects,  the  experimental  findings  of  Morrison  and  Kronauer1,  Morrison  ot. 
al2.  and  Bullock  et.  al!  in  which  the  detailed  structure  of  the  streamwise  component 
of  turbulence  was  investigated  in  pipe  flow.  In  these  experimental  investigations,  evi¬ 
dence  is  presented  for  the  exis'ence  of  a  similarity  family  of  waves  which  extend  from 
the  sublayer  out  into  the  logarithmic  region.  Also  in  that  work,  a  spectral  decomposi¬ 
tion  of  the  velocity  field  in  the  sub-layer  revealed  an  elongated  structure  whose  spamvise 
length  is  close  to  the  generally  accepted  value  for  the  streak  spacing  observed  in  flow 
visualization  studies.  In  addition,  wave  speeds  in  the  sublayer  were  found  to  be  several 
times  higher  than  the  local  mean  velocity. 

These  experimental  findings  have  shown  that  a  horizontal  spectral  decomposition 
of  turbulence  is  remarkably  revealing  about  the  nature  of  turbulent  structure.  It  should 
be  noted,  however,  that  in  the  experiments  cited,  only  the  streamwise  component  of 
turbulence  was  measured.  Strictly  speaking,  however,  to  test  the  Morrison- Kronauer 
wave  similarity  hypothesis,  velocity  measurements  must  be  made  in  the  wave  system  of 
coordinate's.  However,  as  Morrison  and  Kronauer  point  out.  it  is  only  the  wall-normal 
component,  of  turbulence'  which  remains  invariant  in  wave  e'eiordinate's  se>  that  a  strie't 
te'st  of  the  similarity  hypothesis  can  only  be'  made'  by  ele'terminiug  the  structure  of 
this  componc'nt  of  the-  vedeicity.  Experimental  eliflieult ie-s  make'  the'  ele-te-rminatieni  of 
the'  wall-normal  component  very  eliflieult  near  the'  wall.  This  is  e>ne'  motivation  for 
th<‘  ele'te'rminatiem  e>f  the'  de-tails  e>f  the  spectral  strue  ture-  e>f  turbulence-  using  a  elire'et 
nume’rie'al  simulation. 

This  we>rk  is  further  meet ivate’el  by  a  ne'e-.l  for  mem'  ele'taileel  information  em  the- 
nature'  e>f  the'  seeure'e's  of  wall  pre-sslire'  flue't  uat  ions.  Ill  1 1  Ml  ielle'1'.  e't.  al.  1  it  was  shown 
that  a  relafive-ly  low  re'solution  calculatiem  e'.\l:ibits  main-  of  the  feature's  of  the  wall 
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pressure  field  observed  in  experiments.  However,  no  details  of  the  turbulent  velocity 
structure  which  gives  rise  to  the  wall  pressure  structure  were  given.  It  is  well  known, 
however,  that  the  wall  pressure  is  the  result  of  turbulence  interactions  throughout  the 
entire  flow  domain.  Specifically,  the  pressure.  ;>(/,  /).  for  an  incompressible  flow  is  given 
by  the  Poisson  equation: 

V-p(xJ)  =  —  V-(V-T).  (1) 

!> 

where  p  is  the  density  and  T.  the  Reynolds  stress  tensor,  is  given  by: 

T  =  (2) 

and  7  =  i/„f  „  is  the  velocity.  (In  the  current  work,  the  the  subscripts  1.2,  and  3 
will  designate',  for  all  variables,  the  strejunwise.  wall-normal,  and  spanwise  components 
respectively.)  If  the  velocity  is  then  decomposed  into  its  ensemble  averaged  component 
denoted  by  <>.  which  for  a  statistically  stationary  How  is  independent  of  time,  and  its 
fluctuating  components,  we  obtain  for  the  pressure  fluctuation,  ]>'{x.t). 
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where  primes  indicate'  a  turbulent  or  fluctuating  component  and  U\  is  the  mean  ve¬ 
locity  which  is  assumed  to  vary  only  in  the  x->  direction.  It  has  been  argued  by  souk- 
(Kraiclman".  Meecham  and  Tavis1')  that  the'  second  term  on  the  right  hand  side  of 
equat  ion  3.  the  so  called  turbulence  -  t  urbuleuce  (T  —  T)  source',  is  e  >f  see-one  la  ry  impen- 
taiie-e' fee  the' first  .  t he-  t urbuleuce'- -me 'tin  she’ar  ( T  —  .17 )  term.  Chase-'  estimates  that  the 
ratio  e>f  7’  .1/  to  7’ -  T  is  about  2 .  b .  It  is  ge-neTallv  ae-e-e-pte-el.  heme-wr.  that  the  T  M 

se>urce-  is  elominant.  though  perhaps  ne>t  thremghemt  the  e-ntire-  wavenumber  elouiam. 
Sill'-e  the  T  -  M  te-nil  is  <  le  t  ermine  -el  emly  by  t  lie'  Wall- 1  ie  illlial  e'e  till  penie ’lit  e>f  t  ill  bulene'e. 


its  structure  can  be  readily  computed.  This  motivates,  in  part,  the  direction  of  the 
present  work. 

2.  COMPUTATIONAL  PARAMETERS  AND  DATA  ANALYSIS 

The  Navier-Stokes  equations  an'  solved  in  rotation  form  in  a  channel  flow  geometry 
using  a  pseud*  spectral  procedure  which  is  described  in  Orszag  and  Kells8.  In  units  of 
the  channel  half-width,  h  .  the  channel  dimensions  are  5  x  5  x  2  in  the  xj,  X3,  and 
x2  directions  respectively.  The  flow  is  computed  at  G4  x  G4  evenly  spaced  grid  points 
in  the  horizontal  (,i\  —  .r.j )  plain'  and  Go  points,  with  Chebyshev  scaling,  are  used  in 
the  j'9  direction.  The  initial  mean  velocity  profile  is  parabolic  to  which  two  and  three 
dimensional  finite  amplitude  modes  have  been  added.  The  Reynolds  number  (7?)  is 
5000  based  on  h  and  U 0.  the  initial  laminar  centerline  velocity,  and  since  h  and  Uo  are 
defined  as  unity,  R  is  1/;/,  where  r  is  the  kinematic-viscosity.  The  Reynolds  number 
R  is  constant  for  the  duration  of  the  calculation.  The  Reynolds  number  R*  is  defined 
as  u*h/u  where  i /*  is  the  friction  velocity  defined  by  \Jt / p,  and  r  is  the  mean  shear 
stress  at  the  wall. 

As  discussed  in  detail  in  Orszag  and  Patera'1,  the  flow  will  undergo  transition  to 
turbulence  on  a  convective  time  scale.  It  will  then  remain  in  a  quasi-stationary  state 
and  then  decay  to  some  final  state  on  a  viscous  time  scale.  The  final  state  is  determined 
by  the  initial  driving  force  which  is  fixed  in  this  calculation  at  2/7?.  In  the  present  rase, 
a  quasi-stationary  state  is  reached  at  about  th/U (t  —  39.  At  this  time,  data  is  stored 
at  512  time  steps  with  a  sampling  time.  Alh/L'lt,  of  0.0375.  The  computational  time 
step  duration  is  0.00G25.  During  this  quasi-stationary  period.  7?*  varies  from  200  to 
1SG  with  ;m  average  value  of  about  193. 

The  prim  ipal  method  used  to  interpiet  the  data  generated  in  the'  presenr  calcu¬ 
lation  is  spectral  analysis.  The  reader  is  refei  red  to  [2]  for  a  thorough  discussion  of 
th*'  usefulness  of  employing  a  l  ouriei  decomposition  as  opposed  to.  say.  the  spare-tune 
correlation  in  facilitating  the  intei  pi  elation  of  turbulent  structure.  In  the  present  work. 


it  is  natural  to  take  Fourier  transforms  of  the  raw  space-time  data  in  the  homogeneous 
(xi  —  £3)  plane  at  selected  distances,  x2,  from  the  wall.  The  spectral  function  of  interest 
in  the  current  work,  (h, j,  is  defined  by: 


$,.(£, «.\x2,.r'2)  = 


LiL:lT 


<  Ui(k,U),X2)Uj(k,  U>,  X2')  >, 


and. 


(4a) 


—  oc  —  oc  —  oc 

where  i',(  ,r2  )  is  the  Fourier  transform  of  Uj{.v,t ),  A'  is  the  wavenumber  vector  in  the 

horizontal  plane,  a.’  is  the  radian  frequency,  L  is  the  length  (or  width  of  the  domain), 
T  is  the  run  time  for  one  realization,  and  U  is  the  complex  conjugate  of  U.  Here, 
the  overbar  is  defined  by  equation  Oa  with  an  averaging  time.  Ta ,  equal  to  9.C  and  11 , 
represents  the  velocity  with  its  mean,  t/(,  removed. 

In  the  present  work  the  ensemble  average  is  taken  over  two  realizations  each  of 
duration  th/U  —  9.0.  Although  this  gives  a  poor  estimate  of  the  projection  spectra 
defined  below  have-  much  smaller  random  errors.  We  note  that  the  data  set  is  limited 
in  terms  of  total  run  time  compared  to  experiments  so  that  if  a  shorter  averaging  time 
was  chosen  to  increase  the  number  of  realizations,  the  low  frequencies  could  not  be 
resolved.  A  more  detailed  discussion  of  the  limitations  of  the  calculation  is  given  in 

?i3/». 

It  is  also  convenient  to  non-dimensionalize  ‘Iq,  using  the  viscous  length.  /*  —  u/u* . 
viscous  time,  f  —  and  square  of  the  variance.  u,iir 
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where  k*  —  A-/*,  and  u>*  =  u it* .  The  projection  spectra  are  then  defined,  as  in  [2],  by  : 


+  oo 

=  J  tijdkl 


+  OC 

=  J  *ijdk l 


In  this  work,  it  is  most  useful,  for  reasons  given  in  [l]  and  [2].  to  compute  the  spectra 
A ,j  defined  by: 


.\,j(k*.u.%*)  A  *u.-*vl»,y(  A  *.u.-* ). 


a;/(A-;.a;)-a*a^17(A7.a*). 


In  these  expressions,  the  dependence  on  .id  =  ./•■_> //*  and  .r*/  is  implied. 
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3.  RESULTS 

A.  Stationarity  of  the  Flow  and  Turbulence  Profiles 

The  scaling  arguments  given  in  [9]  show  that  the  flow  being  simulated  evolves  on  a 
viscous  time,  1/ v.  That  is,  a  truly  steady  state  cannot  be  achieved  unless  the  integration 
time  is  very  long.  However,  during  transition  to  turbulence  there  docs  exist  a  quasi- 
equilibrium  between  the  inner  layer  and  the  core  of  the  flow.  Classical  arguments  show 
that  under  these  conditions  a  logarithmic  layer  should  exist  simultaneously  with  a  well 
developed  inner  layer.  Though  the  inner  layer  develops  rapidly,  the  core  continues  to 
evolve  slowly.  It  is  during  this  quasi-equilibrium  state  that  the  following  results  have 
been  obtained. 

In  figure  1  (a-d)  we  present  the  mean  velocity  profile,  and  the  streamwise,  nor¬ 
mal,  and  Reynolds  stress  intensity  profiles.  For  each  plot,  there  arc  four  curves  which 
are  obtained  by  performing,  for  a  fixed  distance  above  the  wall,  an  average  over  the 
horizontal  plane  and  a  temporal  average  over  one  quarter  of  the  total  time  used  in  the 
analysis.  This  averaging  procedure  is  represented  by: 


"U'i.M  =  TjT  I 
1 »  •/((> 


<  n(.v,  t, ) 


I  \  r  Li 


<  ii(-r.t)  = 


uu 


where  the  averaging  time.J’,,.  is  4.8  and  each  curve  is  distinguished  by  which  has 
values  39.2.  44.0.  48. S  and  53. G  as  indicated  in  figure  1.  Also,  the  data  in  each  time 
interval  is  sealed  by  the  average  u*  for  that  time  interval. 

The  mean  velocity  profiles  (fig.  1-a)  show  a  well  developed  sublayer  in  the  region 
0.23  <  .)  j  <  3.7.  A  logarithmic  region  exists  for  G5  <  .r.j  <  150.  The'  flow  in  this  region 
is  only  quasi-stationary  as  indicated  by  the*  slow  decrease  in  the  mean  velocity  with 
time.  However,  to  a  good  approximation,  the  expression  n  i/e*  =  .415~ 1  In.r  .J  4-  G.S, 
describes  the  mean  velocity  profile  in  this  region. 
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The  streamwise  and  wall-normal  turbulence  intensity  profiles  are  shown  in  figures 
1-b  and  1-c.  These  results  may  be  compared  with  those  of  Kreplin  and  Eckclmann10 
who  made  measurements  in  a  channel  with  a  Rejuolds  number  based  on  half-width 
and  centerline  velocity  of  3S50  compared  with  the  somewhat,  lower  value  of  3GS0  in  the 
present  study.  The  comparison  is  good  for  i/j "'*/»*  for  <  5  and  for  / a*  for 
x*  <  40.  The  root-mean-scpiare  value  of  a  fluctuating  quantity, /rms,  is  defined  here  as 
frm*  —  s/p.  The  experimental  results  indicate  a  peak  value  for  j\i*  of  2.85  at 
.r*  ~  13  and  a  peak  for  u!')",*/u'r  of  1.0  at  x*,  =  40.  The  calculations  show  peak  values 
for  u\",H /u*  ranging  from  3.04  to  2.2S  at  x*,  ~  Go  and  peak  values  for  u'2"'s /u*  ranging 
from  1.91  to  1.2S.  which  occurs  at  or  near  the  centerline.  In  figure  1-d,  the  U1H2/11*2 
profiles  an'  in  good  agreement  in  the  core  region  with  the  asymptotic  (/?*  ->  00)  limit 
for  the  second  half  of  the  calculation  but  the  stresses  are  too  high  early  in  the  run. 

These  disagreements  between  experiment  and  calculation  are  due  primarily  to  two 
factors.  The  first  is  the  relatively  short  duration  of  the  quasi-stationary  turbulence. 
During  this  time  the  outer  regions  of  the  flow  are  not  well  developed.  This  problem 
can  only  be  remedied  by  computing  much  further  out  in  time.  This  has  been  done 
in  a  recent  higher  resolution  direct  calculation1*  which  also  starts  with  laminar  initial 
conditions.  In  that  calculation,  a  constant,  driving  force  is  chosen  so  that  the  steady 
sta.t.e  wall  shear  stress  achieves  a  value  significantly  higher  than  the  initial  laminar 
value.  At  the  end  of  the  calculation  the  peak  in  u™,a/u*  occurs  at  x*  ~  23.  The 
run  time  of  the  higher  resolution  calculation  will  be  significantly  increased  to  see  if 
convergence  with  experiment  is  obtained.  We  also  note  that  in  [11]  it  is  shown  that  the 
incorporation  of  subgrid  scale  models  obtained  from  a  renormalization  group  analysis 
gives  turbulence  peaks  much  closer  to  the  wall  in  a  relatively  short  time  starting  from 
laminar  initial  conditions. 

The  second  source  of  disagreement  is  the  lack  of  good  spanwise  resolution.  Orszag. 
et.  al.  have  shown  that  good  spanwise  resolution  is  crucial  and  this  is  due  to  the 
necessity  of  capturing  the  near  wall  dynamics  of  the  bursting  process.  Recent  results 


of  Leighton12,  however,  using  the  identical  data  set  used  in  obtaining  the  results  of  this 
paper,  show  that  the  bursting  process  is  at  least  nominally  captured  by  the  current 
spamvise  resolution  of  about  13  viscous  units.  With  these  qualifications  in  mind,  the 
structure  of  the  turbulent,  field  will  now  be  described. 

B.  Visual  Observations  and  Spectra 

Certain  structural  features  of  the  streamwise  and  wall  normal  velocity  fields  can  be 
determined,  at  least  qualitatively,  by  viewing  contour  plots  of  the  velocity  in  horizontal 
planes  at  a  fixed  instant  in  time.  In  figures  2(n-d)  and  3(a-d)  we  have  exhibited  such 
plots  for  the  u  i  and  t/2  fields,  respectively.  For  the  i/j  component,  we  plot  the  contours 
( t/ 1  —  <  </j  >/, )  —  u\"‘s  and  (  u  i  —  <  u  i  >/,)  —  (  0.25  x  u  i'"'*  )  and  for  </>  w<'  plot  the  contours 
u  >  +  i/ and  u  >  +  (0.25)  X  That  is.  the  streamwise  velocity  contours  correspond 

to  fluctuations  in  the  opposite  direction  to  the  mean  flow  and  the  wall-normal  velocity 
contours  represent  fluctuations  away  from  the  wall.  These  contour  levels  are  chosen  to 
correspond  to  low  speed  regions  of  the  How  which,  near  the  wall,  are  normally  observed 
in  How  visualization  studies  as  streaks.  In  effect,  two  such  plots,  one  of  t/  i  and  one  of  t/2 
at  same  elevation  above  the  wall,  will  give  instantaneous  contours  of  second  quadrant 
Reynolds  stresses. 

In  Hgure  2-a.  (.r*,  =  2.1 1 )  we  note  that  theiq  velocity  structure  is  clearly  elongated 
in  the  How  direction.  The  corresponding  plot  of  u  <  in  figure  3-a  also  shows  an  elongated 
structure  and  is  very  strongly  correlated  with  the  corresponding  </|  structure.  However, 
the  normal  velocity  structure  is  distinctly  different  from  the  streamwise  structure  in 
that  the  u  •  contouis  are  not  strongly  elongated  in  the  How  direction  and  appear  to  be 
distinctly  finer  grained.  At  =  27.7.  we  observe  a  u  \  structure  at  a  spamvise  location 
of  about  3.5  that  extends  the  entire  length  of  th<’  How  domain.  Xo  such  structure  is 
clearly  evident  in  the  n  ,  field.  At  ./A  =  9-1.75.  the  edge  of  the  logarithmic  layer,  only 
ii  |  st  met  tires  of  spa  n  wise  wavelengt  h  on  t  he  order  of  the  box  dimensions  exist  whereas 
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ill  the  corresponding  u  >  plot  the  structures  are  not  nearly  as  large.  We  also  note  a  lack 
of  any  significant  correlation  between  u  i  and  u  ,  at  this  level  above  the  wall. 

Before  discussing  the  spectral  structure  of  these  fields,  it  is  appropriate  at  this 
point  to  indicate  the  limitations  of  the  present  calculations,  in  addition  to  the  problems 
associated  with  non-stationarity  and  the  early  state  of  development  of  the  turbulence 
discussed  earlier.  These  constraints  are  those'  imposed  by  the'  computational  domain 
size-  and  grid  spacing,  which  determine  the  lowest  and  highest  resolvable'  wavenumbers, 
and  total  computational  time'  and  time  step  size,  which  eletermine  the'  lowest  and  highest 
resolvable'  fre'que'iivies.  In  the'  e-urre-nt  cale-ulafiem.  the1  griel  re'solutiou  in  'he  horizontal 
plane1  is  about  15.2  in  vise-tms  units  anel  the'  heuizemtal  plane'  is  abemt  975  x  975  in 
the'se'  same'  units.  These'  elimensiems  give-  0.00G4  anel  0.205  .  fe >r  the'  lowe'st  anel  highe’st. 
wave'immbe'rs.  }>■*,  respe'ctivedy.  Similarly,  the-  total  run  time'  in  vise’ems  units  is  14G  anel 
the'  sampling  time  is  0.2S5.  This  give's  0.0-43  and  11.0  for  the-  lenve-sf  anel  highest  raelian 
fre'ejue'iicie'S,  re-spe-e-tive-ly. 

Henvewe-r.  the'se  limits  em  the'  range  of  the-  data  in  wave'iiumbe'r- frequency  ceiordi- 
nates  were*  moelifiee!  because'  of  computer  memory  limitations  in  the1  preice'ss  of  cemi- 
puting  spe-e  t l ;i  freim  the  raw  elata.  These  limits,  again  in  viscous  units,  are  .00044  anel 
0.0900  for  the'  lowe-st  ;mel  highe'st  stre-amwise'  ;mel  spanwise-  wave-numbers  anel  0.0SG1 
anel  11.0  for  the-  lowe-st  anel  highe'st  raelian  fre-ejue-neies.  The'  analysis  is  the-re-feue'  cem- 
fineel  tej  approximate-ly  ;i  one  orde-r  e >f  magnituele-  variatiem  in  wavi'niunbe-r  anel  twei 
oreh-rs  of  magnituele-  in  fre-eiue-ne-y.  We-  also  ne>te'  that  spanwise-  strue-ture-s  with  wave- 


le-ngths.  A*.  e>n  the-  e>rele-r  of  100.  the-  so  e-alle-et  sublaye-r  stre'ak  spacing,  a l e-  e-xpe-cteel 
to  appe  al  at  A-.*  -----  0.0G3  see  that  sue’h  information  is  within  (though  aehnitte-elly  at  the- 
very  high  enel)  of  the-  e'ale’ula t e-el  spectra.  We-  also  eibse-rve-  that  in  e-xpe-rime-nts' -■■  .  at 
least  a  two  eneler  e it  magnitude  range-  was  cove-re-el  in  /.'*  vC*  e-oorelinate's.  4  his  allows  a 
re-solution  of  not  only  the-  small  structures  near  the-  wall  but  large-  strue-ture-s  eibservexl 
in  the-  leigaritlnnie-  portiem  of  tl-e  velocity  preifile-.  In  this  region,  significant  strue-ture' 
is  obse-rve-el  for  raelian  irei|uencii's  lie-low  10'"".  It  is  elearly  not  possible  to  investigate 


4 


i  t  *  a-4  •  I  I  )  '1.1  l.t'l.t  i.i'4,1  tit  tat  »-*  |.l  I  I  |.t  *-  4  a .  ■  »  t'ta*  ■  -  >  » _ I '  1.1  §-t  O  a. I  1 .  > 


^a'J.  -la'  J.'  t.  • 


this  region  of  the  spectrum  in  the  present  study.  However,  certain  qualitative  behavior 
of  the  spectral  structure  can  be  discerned. 

With  these  limitations  in  mind,  examples  of  the  projection  spectra  .  An  and  Ajj. 
are  shown  in  figures  4  ( a  - 1 )  and  5  (a-d)  respectively.  In  these  plots  the  spectra  displayed 
are  those  for  which  ,r>*  =  .r >*'.  In  the  sublayer,  at  .r*  =  5.S4.  the  t/j  field  exhibits 
a  definite  peak  in  the  k*  —  u.1*  spectrum,  but  the  convective  ridge  is  not  well  defined. 
The  k\  —  A-3  spectrum  exhibits  a  well  defined  peak  at  />•■*  ~  0.09  and  k'  ~  0.02G.  These 
values  correspond  to  a  spanwise  wavelength  of  09.7  and  a  streamwise  wavelength  of 
243,  indicating  structures  which  are  clearly  elongated.  It  would  be  tempting  to  state 
at  this  point  that  this  peak  in  spectrum  corresponds  to  the  streaky  structure  observed 
experimentally  by  many  investigators.  The  spanwise  wavelength  of  09.7.  however,  is 
considerably  smaller  than  the  generally  accepted  value  of  about  100.  Morrison  and 
Kronauer"  obtained  Aj  =  130  in  their  experiment  using  the  same  spectral  analysis.  On 
the  other  hand,  some  investigatiors  have  found  values  for  A3  less  than  100,  notably  the 
measurement  of  Runstadler1 '.  Though  we  cannot  draw  a  definite  conclusion  as  to  the 
origin  of  this  peak  at  the  present  time,  it  is  encouraging  that  the  k^  —  k%  spectrum 
corresponds  to  the  visual  evidence  of  an  elongated  structure  shown  in  the  contour  plots. 
More  detailed  analysis  of  higher  resolution  data  is  currently  in  progress  to  determine 
more  clearly  the  origin  of  this  peak. 

At  .1  •*  =  27.7  we  set'  that  the  /.•J'— W  spectrum  of  the  streamwise  velocity  component 
has  a  clearly  defined  convective  ridge.  It  is  also  clear  from  the  k*  —  k‘f  spectrum,  that 
the  fluctuation  energy  begins  to  shift  to  lower  wavenumbers  at  this  distance  from  the 
wall.  At  ,r*2  —  04.  the  bottom  edge  of  the  logarithmic  layer,  this  shift  is  evident.  A 
spectral  peak  is  clearly  defined  in  the  /c j  —  W  spectrum  at  /.*  =  0.02S  or  A3  =  223. 
There1  is  a  further  spectral  shift  at  .1*,  =  100.4.  the1  outer  part  of  the1  logarithmic  layer, 
to  about  /e*  =  0.019  or  A3  =  331.  This  wavelength  corresponds  roughly  to  that  of  the 
width  of  the  channel.  2//.  A  transition  from  length  scales  on  the  order  of  the  streak 
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spacing  to  that  of  the  channel  width  seems  to  occur  smoothly  as  we  move  across  the 
buffer  layer. 

This  transition  of  length  settles  is  also  clearly  evident  in  the  transition  of  the 
local  convection  velocities  as  shown  in  figure  G.  In  this  figure  we  have  computed  the 
convection  velocitv  of  the  sfreamwise  fluctuations  at  several  locations  above  the  wall  by 
selecting,  at  various  wavenumbers,  the  frequency  which  corresponds  to  the  maximum 
in  the  k *  -w*  spectrum.  An  estimate  of  the  convection  velocity  can  then  he  made  using 
the  average  value' of  c*  =  k^ /^'*  at  several  locations  along  the  convective  ridge.  In  figure 
G.  wo  plot  r* .  (the  convection  velocity,  c.  normalized  with  a*  )  ami  the  ratio  cf  <  1/1  > 
where  <  u  \  >  is  the  local  st teamwise  velocity  which  is  obtained  by  averaging  both 
horizontally  and  over  the  total  time  of  the  run.  The  local  convection  velocity  is  clearly 
much  larger  than  the  local  mean  speed  in  the  sublayer.  Tins  is  not  inconsistent  with 
experimental  results  [2]  which  indicate  that  the  convection  velocities  normalized  with 
the  local  mean  velocity  vary  from  about  5  at  ,rj  =  1.52  to  2  at  ,rj  -  5.93.  The  values 
obtained  here  tire  somewhat  higher,  ranging  from  G.S  at  r*,  —  2.1  to  2.7  tit  .v',  —  5.04. 
The  calculations  also  show  that  the  convection  velocities  in  the  sublayer  tire  remarkably 
constant  tit  about  c*  =  14.  This  is  higher  than  the  experimental  result  of  [2]  in  which 
c*  is  found  to  be  about  S.O.  but  the  experiments  also  show  that  this  speed  is  unchanged 
throughout  tin' sublayer.  For  .id,  >  10.  the  calculation  shows  that  r*  gradually  increases 
to  about  19  in  the  logarithmic  layer  which  is  only  slightly  above  the  local  moan  flow 
velocity.  All  of  these  results  are  consistent,  at  least  qualitatively,  with  the  experimental 
results  though  it  is  not  possible  on  the  basis  of  the  current  computational  evidence  to 
confirm  or  reject  the  hypothesis  proposed  m  [2 j  that  such  high  speeds  in  the  sublayer 
are  in  fact  due  to  a  wave  phenomenon. 

He]  >  resent  at  i  ve  spectra  for  the  u  <  fluctuations  shown  in  figure  5  are  dillerent  in 
some  respects  fj om  the  n  \  sped i  a.  I  his  should  cei  t a mly  be  expected  from  the  obvious 
di  defences  shown  m  the  coi  it  om  plots.  I  lie  most  not  able  <  h  defence  is  the  absence  of  an 
energy  peak  at  wavelengf  lis  cm  i  esj >,  >m line,  to  t  he  sublayer  st leak  spacing.  I  he  -  k\ 
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spectrum  for  x*2  =  11.4  indicates  a  possible  peak  at  a  coordinate  much  higher  than 
that  found  the  uj  field.  Also,  at  x%  -  G4,  the  edge  of  the  logarithmic  layer,  the  u> 
energy  begins  to  shift  noticeably  toward  large  length  scales.  The  energy  has  clearly 
shifted  to  scales  on  order  of  the  channel  size  at  .v*,  =  1GG.3.  It  therefore  appears  that 
the  u 2  field  has  generally  a  smaller  structure  than  u  \  and  this  is  completely  consistent 
with  visual  observation.  The  author,  however,  is  not  aware  of  any  experimental  data 
on  the  wall  normal  velocity  which  can  be  used  for  comparison. 


As  discussed  in  [15]  and  [1G],  we  may  expect,  from  theoretical  considerations  . 
certain  differences  between  the  structure  of  the  wall  pressure  field  in  a  channel  flow 
and  its  structure  in  a  boundary  layer  flow.  In  particular,  in  strictly  incompressible 
How  we  expect  the  spectral  density  of  the  pressure  field  to  vanish  as  the  wavenumber 
approaches  zero.  This  dot's  not  appear  to  Ire  the  case  in  a  channel  How.  Nevertheless, 
we  should  not.  expect  significant  differences  in  the  overall  structure  of  the  wall  pressure 
field. 

In  figure  7  (a-c),  we  show  the  space-time  correlations  of  the  pressure  field  obtained 
in  the  present,  calcultion.  These  correlations  are  defined  by: 
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In  these  expressions,  the  spectrum  of  the  wall  pressure,  vI'/)/),  is  defined  in  an 
analogous  fashion  to  the  spectrum  of  the  velocity  field.  'I'(J.  given  in  equation  7.  It 
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is  convenient  to  examine  the  pressure  field  in  space-time  coordinates  (rather  than  in 
transform  space)  for  present  purposes  simply  because  in  these  coordinates  the  gross 
behavior  of  the  results  are  more  easily  interpreted  and  compared  with  the  available 
experimental  results.  (In  [1G]  the  pressure  field  structure  is  given  in  spectral  coordinates 
and  the  ratio  of  the  root  mean  square  pressure  to  the  wall  shear  stress  is  found  to 
be  3.915.)  The  function  Rjtl,(  A.r* ,  At* )  shows  a  clearly  defined  convective  ridge.  A 
line  showing  the  approximate  position  of  this  ridge  shows  clearly  that  the  convection 
velocity,  c*  =  dAx* /dAt* .  (that  is  the  slope  of  the  line  defining  this  ridge),  is  not 
constant.  The  convection  velocity  at  very  small  streamwise  spatial  separations  is  about 
12.1  and  at  very  large  separations  (say  A.rJ  >  300)  it  is  about  13.4.  The  ratios  of  these 
velocities  to  the  centerline  velocity  are  0.G4  at  small  separations  and  0.70  at  large 
separations.  These  results  arc  certainly  in  qualitative  agreement  with  the  results  of 
Blake1'.  Dull 1 8 ,  and  Willmarth19.  It  is  also  consistent  with  the  idea  that  since  the 
pressure  field  is  driven  by  sources  that  exist  throughout  the  flow  domain,  one  expects 
smaller  convection  speeds  at  small  separations  (say.  separations  that  are  at  least  an 
order  of  magnitude  smaller  than  the  largest  scales  of  the  How)  since  this  region  of  the 
space-time  correlation  may  be  more  representative  of  structures  whose  centers  are  close 
to  the  wall  where  local  convection  speeds  art'  small. 

Also,  from  elementary  theoretical  considerations  we  must  expect  that  if  the 
turbulence-mean  shear  source  term  is  dominant,  that  the  factor  A-i/|A'|  =  cos#.  where 
k  =  | A- 1  and  9  is  t lit'  angle  between  k  and  the  How  direction,  will  dominate  the  spec¬ 
tral  characteristics.  This  struct  me  is  indeed  observed  in  the  experiments  of  Bull1  ‘  who 
measured  /?/)/(( A.r*.  A.r* ).  For  a  constant  correlation  value  of  0.1  he  obtained  a  contour 
whose  elongation  in  the  spanwise  direction  was  about  2.7  times  that  in  the  streamwise 
direction.  The  present  calculations  give  a  ratio  of  about  3.3.  In  addition,  wo  note  that 
the  calculations  give  approximately  concentric  correlation  contours  for  /?;i/,(  A.r* .  A.r* ) 
at  small  separations  indicating  a  nearly  isotropic  small  scale  field.  This  is  also  in  excel¬ 
lent  agreement  with  experimental  results.  Finally,  the  7?;1/1(  A.r*.  Af  * )  contours  show  no 
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indication  of  spanwise  convection  and  are  almost  perfectly  symmetric  about  At*  =  0. 
In  general,  the  qualitative  agreement  between  these  wall-pressure  results  and  available 
experiments  is  excellent,  but  the  nature  of  the  sources  of  the  pressure  field  can  only  be 
obtained  by  computing  the  dependence  of  the  velocity  correlations  on  the  wall-normal 
coordinate.  To  this  end.  preliminary  results  on  the  wall-normal  structure  of  i/2  have 
been  computed  and  are  presented  next. 

D.  Wall-normal  Structure  of  the  Vertical  Component  of  Velocity 

As  indicated  in  the  introduction,  the  structure  of  the  uo  velocity  component  de¬ 
termines  the  structure  of  the  T-M  source  of  the  wall  pressure  field.  As  a  preliminary 
indication  of  this  structure,  we  compute  both  the  broadband  correlation  coefficient, 
.iy>.  and  the  correlation  coefficient  in  narrow  frequency  bands,  n 22.  These  functions 
are  defined  bv: 
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where  ‘I’jj  is  defined  by  equation  (  1 ).  \\  < •  also  compute  the  phase  difference,  tf. 


given  bv: 
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In  figure  (Sa)  we  plot  .if,  )■  For  each  curve,  is  fixed  at  11.4.  C4. 03. 129.31. 

and  1G6.39  in  viscous  units  and  .i  >  i-s  varied.  The  asymmetry  of  ,i>2  is  clear.  That  is,  for 
each  curve,  the  rate  of  decay  of  the  correlation  coefficient  is  larger  for  r>  <  ,i  '2 .  We  also 
note  that  is  positive  across  the  entire  width  of  the  channel.  These  two  characteristics 
arc  qualitatively  consistent  with  the  experimental  measurements  of  Comte- Bellot."0.  For 
clarity  at  smaller  separations,  this  same  data  is  shown  on  a  semi-log  plot  in  figure  S(b). 
In  figure  8(c),  the  same  data  is  is  plotted  against  .v> /.r”  where  x’2  ^  is  the  fixed 
reference  coordinate  for  each  plot.  It  is  clear  that  the  broadband  correlations  cannot 
be  forced  to  collapse  over  the  entire  width  of  the  channel  although  there  appears  to  be 
a  limited  region  of  collapse  for  .v2*  =  G4. 03. 129.31  and  1GG.39  for  <  1. 

The  wave  similarity  hypothesis  proposed  in  [1]  and  [2],  however,  suggests  that  u; 
may  be  list'd  to  constrain  |  />•  |,  at  least  approximately.  This  led  Bullock,  et.  al . [3]  to 
propose  that  02  >  could  be  reduced  to  a  dependence  on  only  two  parameters, 
and  .  With  this  as  motivation  we  plot,  in  figure  9(a-c),  022  for  three  values  of 

u-'j'V  ^ ,  though  the  value  of  this  parameter  could  not  be  held  absolutely  constant  because 
of  the  discrete  nature  of  the  data.  The  collapse  is  not  nearly  as  good  as  the  experimental 
data  given  in  [3]  but  part  of  the  reason  for  this  can  clearly  lie  attributed  to  the  lack  of 
sufficient  data  needed  to  reduce  statistical  scatter  in  the  spectral  estimates.  However, 
it  is  clear  from  these  plots  that  the'  correlation  length  is  decreasing  with  increasing 
v* which  is  consistent  with  the  experimental  results  given  in  [3].  In  figure  [10]  we 
plot  the  phase  change,  ffi  given  by  equation  12,  against  ,r  >  for  .v',‘ ^  =  G4.05  for  four 
different  frequencies.  This  again  shows  a  typical  asymmetry  with  respect  to  the  side 
of  the  channel  above  and  below  the  fixed  reference  point.  We  also  note  larger  phase 
differences  for  larger  frequencies.  These  results  are  also  qualitatively  consistent  with 
the  results  of  [3],  though  in  that  work  only  the  streamwise  component  of  turbulence 
was  measured. 
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4.  SUMMARY  AND  CONCLUSIONS 


The  wavenumber-frequency  structure  of  a  quasistationary  turbulent,  flow  has  been 
computed  from  data  generated  from  a  direct  simulation  in  a  channel.  Although  the 
turbulence  is  in  an  early  stage  of  development.,  certain  features  of  the  velocity  structure 
compare  favorably  with  experimental  results.  In  particular,  the  spanwise  periodicity 
of  the  streamwise  velocity  component  in  the  sublayer  is  clearly  evident  in  contour 
plots  of  the  field  and  in  the  corresponding  k\  —  k:i  spectrum  which  indicates  a  spanwise 
wavelength  of  C9.7  (in  viscous  units)  for  this  structure.  The  local  wave  speeds  computed 
from  the  k\  —  jj  spectrum  of  the  streamwise  velocity  component  are  as  much  as  C.S  times 
the  local  average  streamwise  velocity  in  the  sublayer.  The  ratio  of  the  wave  speed  to 
the  local  mean  speed  decreases  across  the  buffer  layer  and  is  approximately  equal  to  one 
in  the  logarithmic  layer.  Also,  the  ratio  of  the  wave  speed  to  friction  velocity  remains 
constant  at  about  14  throughout  the  sublayer.  In  general,  these  results  are  consistent 
with  the  experimental  results  described  in  [l]  and  [2],  The  spectra  and  contour  plots 
for  the  wall-normal  component-  of  velocity  indicate  a  structure  which  is  more  finely 
grained  than  that  of  the  streamwise  velocity  component.  Additionally,  the  wall normal 
component  lticks  a  clearly  defined  periodic  spanwise  structure  in  the  sublayer. 

The  space-time  correlations  of  the  wall  pressure  field  reveal  features  which  are  re¬ 
markably  close  to  experimental  observations.  First,  wave  speeds  at  small  separations 
are  somewhat  lower  than  speeds  at  large  separations,  which  is  consistent  with  most  ex¬ 
perimental  results.  Secondly,  the  spanwise-streamwise  correlations  reveal  an  isotropic 
structure  at  small  separations  and  highly  anisotropic  structure  at  large  separations. 
At  larger  separations  the  spanwise  length  scale  is  about  3  times  the  streamwise  scale. 
This  is  consistent  with  experiments  as  well  as  with  elementaiy  theoretical  consider¬ 
ations.  Also,  the  wall  normal  structure  of  t/_>.  which  determines  the  structure  of  the 
turbulence-mean  shear  source  of  wall  pressure,  is  computed.  An  asymmetry  exists  in 
the  broadband  correlation  coefficient  with  respect  to  the  fixed  location,  r','  * .  Follow 
mg  the  conjecture  postulated  m  [3],  the  coi relation  coefficients  computed  in  narrow 
frequency  bauds  are  plotted  against  .ty./.i  V  ^  with  ^  fixed.  A  fair  collapse  of  the 
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(lata  is  obtained  with  a  more  rapid  decay  of  the  correlation  for  higher  values  of  u.\r2  . 

The  statistical  scatter  in  these  narrow  hand  correlations  is  considerable. 

Although  several  aspects  of  this  work,  including  a  nonstationary  turbulence  and 
the  limited  range  of  the  results  at  low  frequecies  and  wavenumbers,  limit  its  usefulness 
in  terms  of  comparison  with  the  details  of  experiments,  the  results  are  close  enough  to 
experiments  to  warrant  the  use  of  direct  simulation  data  to  elucidate  the  nature  of  the 
sources  of  the  wall  pressure  field.  In  particular,  data  obtained  from  higher  resolution 
direct  simulations  will  be  used  in  the  future  to  study  the  detailed  structure  of  both  the 
turbulence-  mean  shear  and  turbulence  -  turbulence  sources  of  the  wall  pressure  field 
and  relevant,  pressure- velocity  correlations. 
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Fig.  1  —  (Continued)  Global  turbulence  properties  obtained  from  the  averaging  detined  by  equation 
9.  Curves  are  distinguished  by  the  lower  limits  of  the  time  integration  given  by:  q/„  =  39.2.  .4. 
tn  =  44.0,  +,  r0  =  48. H,  X  ,  f(l  =  53.6. 
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Fig.  1  —  (Continued)  Global  turbulence  properties  obtained  from  the  averaging  defined  b\  equation 
9.  Curves  are  distinguished  by  the  lower  limits  of  the  time  integration  given  by:  =  39.2.  A. 

t()  =  44.0,  +  ,  /„  =  48.8.  X  ,  f„  =  53.6. 
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He  I  (Continued)  Global  turbulence  properties  obtained  from  the  averaging  defined  by  equation 
9  Curves  are  distinguished  by  the  lower  limits  of  the  time  integration  given  by:  =  39.2.  A. 
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(Continued)  Contours  of  constant  wall  normal  velocity.  Contour  levels  given  by 
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(c)  versus  noiuliinensioiuil  wall  distance. 


Fig. 


—  (Continued)  Broadband  correlation  function  for  the  wall-normal  velocity.  Reference  distances 
from  the  wall  in  viscous  units,  — ,  11.4, - .  64.0,  — - ,  129.31, - ,  166.39. 
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